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ABSTRACT 

Prominences and cavities are ubiquitously observed together, but the physical link between these 
disparate structures has not been established. We address this issue by using dynamic emission in the 
extreme ultraviolet (EUV) to probe the connections of these structures. The SDO/AIA observations 
show that the cavity exhibits excessive emission variability compared to the surrounding quiet-sun 
streamer, particularly in the 171A bandpass. We find that this dynamic emission takes the form of 
coherent loop-like brightening structures which emanate from the prominence into the central cavity. 
The geometry of these structures, dubbed prominence horns, generally mimics the curvature of the 
cavity boundary. We use a space-time statistical analysis of two cavities in multiple AIA bandpasses 
to constrain the energetic properties of 45 horns. In general, we find there is a positive correlation 
between the light curves of the horns in the 17lA and 193A bandpasses, however the 193A emission is 
a factor of 5 weaker. There is also a strong correlation between structural changes to the prominence 
as viewed in the He II 304A bandpass and the enhanced 17lA emission. In that bandpass, the 
prominence appears to extend several megameters along the 171 A horn where we observe co-spatial, 
co-temporal 304A and 17lA emission dynamics. We present these observations as evidence of the 
magnetic and energetic connection between the prominence and the cavity. Further modeling work is 
necessary to explain the physical source and consequences of these events, particularly in the context 
of the traditional paradigm: the cavity is under dense because it supplies mass to the over dense 
prominence. 



1. INTRODUCTION 

Prominences are regions of cool, condensed plasma sus- 
pended in the solar corona. Our understanding of promi- 
nence physics has been divided into two distinct focuses: 
the topology of the magnetic support within the promi- 
nence and the mass source for the prominence itself. 
Through the many iterations of experiments, there has 
been an overlooked constraint on the prominence sys- 
tem, coronal cavities. Coronal cavities are the elliptical 
regions of hot (approximately 1-2 MK) plasma which sur- 
round the prominence. Cavities have been observed to be 
density depleted as compared to the surrounding corona. 
Rudimentary calculations, based on cavity volume and 
density, seem to indicate that the mass reduction in the 
cavity is not enough mass to pop ulate the prominence 
( |Saito fc Tandberg-Hanssen||1973i . 

An important aspect of the cavity is the strongly con- 
trasted boundary that is observed between the cavity 
and the surrounding streamer. In the magnetically domi- 
nated corona, strong gradients in plasma properties (that 
would produce gradients in emission) must be related to 
the magnetic field and more specifically to distinct flux 
systems. How is the cavity boundary related to the pres- 
ence of stressed fields? This question is central to our 
approach to studying the prominence-cavity system. 

Much research has been conduct ed attempting to di- 
agnose the properties of the cavity. Gibson et al. (20101 
presents a mult i wavelength morphology study and con- 
tains the most comprehensive reference list. Cavities 



1970 Fuller fc: Gibson|2009 ) and these data measure un- 
equivocally that cavities are a density- depleted feature . 
Studies in X-ray and EUV observations ( Serio et al. 1978 1 
have found that cavities do emit at coronal tempera- 
tures. Later studies in those spectral regi mes uncovered 
that substructure was c ommon in cavities (Hudson et al. 
1999; [Reeves et alT 2012 ), which raises the question: what 
causes the thermodynamics properties of the cavity to 
vary on smaller scales? 

In this article, we detail how dynamic EUV emission 
in the cavity can be used to extract structural and en- 
ergetic information of the prominence-cavity superstruc- 
ture. The layout of the article is as follows: Section 2 
describes generalized observations of dynamics in promi- 
nences and cavities. Section 3 characterizes the emission 
of prominence horns, and Section 4 discusses the conclu- 
sions we are able to draw from these observations as a 
whole. 

1.1. Prominence- Cavity Dynamics 

Prominenc es have long been observed to be dynamic 
phenomena (Rudnick 1934) . The sem inal observations 
are described in (Zirker et al. 19981). Doppler mea- 



were originally seen in white light scattering ( Waldmeier 
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surements are used to show that bulk motions in the 
prominence can be aligned but anti-parallel despite spa- 
tial thread separations of only few hundred kilometers 
(thread widths are on the order of the resolution ele- 
ment « 100 km). This implies a highly coherent mag- 
netic field but with distinctly non-coherent thermody- 
namic processes. This study did not include any coro- 
nal datasets. Higher resolution observations show that 
in addition to thread-moti ons, bubbles are als o a dy- 
namic prominence feature ("Berger et al. 2011). Bub- 
bles are believed to be low density, magnetically- isolated 
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structures which rise through prominence material buoy- 
antly. The prominence material passively shifts around 
the bubbles. Both bubbles and thread motions operate 
on sub-megameter spatial scales, but despite these small 
scale dynamics, the large scale prominence structure re- 
mains intact. Although, both of these observations high- 
light the complex interaction of the magnetic field and 
plasma in the prominence, they both neglect to consider 
the larger role the cavity plays in the dynamics. 

The first observations o f Dopp ler flows in cavities was 
reported in |Schmit et al. ( 2009[ ). The coronal observa- 
tions detailed in that article iaclced the spatial and tem- 
poral resolution of the prominence observations. How- 
ever, it did find large-scale, long-duration flows to be 
present and common in cavities. The source of these 
flows was not determined, but more importantly, the di- 
rect conclusion is that cavities which are observed to sta- 
bly exist for weeks must also be regarded as a structure 
with internal dynamics, analogous to the prominence. 

To extend this work, we have observed cavities in the 
Solar D ynamics Observatory Atmospheric Imaging As- 
sembly ( [Lemen et"al]|2012[ SDO/AIA) dataset. Al- 
though, the AlA dataset lacks spectral information to 
derive Doppler flows, it makes use of high resolution, 
high cadence bandpass imagers which span a broad 
range of temperature sensitivities. With this dataset, 
we can observe the changes in emission which are indica- 
tive of changes in the thermodynamics properties of the 
prominence-cavity plasma. Towards this aim, we are in- 
terested in one particularly relevant paradigm: is a mass 
exchange responsible for moving plasma from the under- 
dense cavity into the over-dense prominence? There are 
various methods to describe this pro cess but the i r com- 
monality lies in thermal instability (Field 1965 Meer- 
son|| 19961). A coronal loop is perturbed such that a seg- 



ment of the loop radiates more energy than the combina- 
tion of heating and conduction input into the segment. 
This causes a runway cooling effect in which the coro- 
nal plasma condenses as it cools to reach a new thermal 
equilibrium in a chromosphere-like condensation state. 
Catastrophic cooling has been related to two diff erent 
processes: shearing field lines (IChoe fc Lee 19921 and 



chromospheric heating (Karpen et al 



2DDS) 



This article focuses on diagnosing the dynamics of the 
prominence-cavity. We establish with this article that 
there is a magnetic and energetic connection between 
the prominence and the cavity. A subsequent article by 
these authors will compare the observed dynamics to the 
hydrodynamic models of catastrophic cooling. This fol- 
lowup article will address the physical source and conse- 
quence of that energetic connection. 

2. EUV DYNAMICS IN THE PROMINENCE CAVITY 

We use the SDO/AIA dataset to observe dynamics in 
the prominence-cavity system in the EUV. SDO/AIA im- 
ages the Sun in eight EUV and UV bandpasses which 
are centered on strong emission lines. The strength of 
these emission lines vary with temperature so the AIA 
datasets present us with a qualitative temperature di- 
agnostic tool. Of the 8 bandpasses on AIA, four of the 
bandpasses are useful for prominence-cavity diagnostics: 
304A, I7lA, I93A, and 21lA. 

The 304A bandpass is centered on the singly-ionized 
Helium Ly-a transition, although it also includes a Si 



XI line which is generally two orders of magnitudes dim- 
mer. He II 304A emission is believed to primarily occur 
between 0.5— 1.0 x lO^K although extended optically thin 
emission throughout the transition region is also possi- 
ble. The 171 A bandpass is centered over a blend of Fe 
IX and Fe X emission which is strongest at a tempera- 
ture of 1 MK. These are the strongest lines by an order 
of magnitude but there is a weak transition region com- 
ponent between 2.5-5x10^ K formed by O V-IV. The 
I93A bandpass centers on an Fe XII transition which is 
characteristic of 1.5 MK plasma. The 21lA bandpass is 
centered on a Fe XIV transition at 211. 3A. In ionization 
equilibrium, Fe XIV reaches a population maximum at 
2 MK. In terms of spectral blends, the 21lA bandpass 
is fairly pure, but its diagnostic value is damped by its 
signal strength. 

SDO/AIA allows us to image the solar plasmas with 
the temperature range of 0.05-2 MK. In practice, the ob- 
servations provide a convolution of many different tem- 
perature structures which complicates our interpretation. 
Moreover, the data contains instrumental effects and ob- 
servational uncertainties, which we must derive and un- 
derstand to accurately state our measurements. 

Throughout this section, we will be discussing two par- 
ticular prominence-cavity datasets. The first is a 30- 
hour dataset collected between 28-Feb-2011 09:00UT and 
l-Mar-2011 15:00UT, hereby known as Dl. The cav- 
ity of interest is off the northwest limb at a latitude of 
60° . Based on the duration of visibility, the longitudinal 
length is estimated at 52° or 320 Mm. The second is a 
36- hour dataset collected between lO-Aug-2011 00:OOUT 
and ll-Aug-2011 12:00UT, hereby known as D2. This 
cavity is off the southeast limb at a latitude of -51°. The 
longitudinal length is estimated at 38° or 340 Mm. Each 
dataset is composed of full resolution data at 6 minute ca- 
dence. Each image in the dataset is the result of median 
filtering 5 sequential images to increase signal-to-noise 
ratio. 

The first steps to understanding the role of dynamics 
in cavities is to understand statistically how the cavity 
differs from the surrounding streamer. 

2.1. Characterization of Prominence- Cavity Structure 

in EUV 

Figurefl] displays the Dl cavity in five AIA bandpasses. 
The 304A bandpass displays the cool, dense prominence. 
This particular prominence extends 18 Mm above the 
limb and has an apparent width of 30 Mm. There is no 
clear separation between the limb and lower boundary of 
the prominence. 

Figure lb shows the prominence in 13lA. The 
I31A bandpass data theoretically contains emission from 
a FeVIII line (blended with the primary Fe XXIII line). 
The I31A bandpass simply does not have enough signal 
to aid in this analysis. 

The signal in the 17lA bandpass is much stronger than 
I31A. The prominence appears in emission. This emis- 
sion is peculiar but ubiquitous in prominences. The peak 
ionization temperature for Fe IX is 8 times higher than 
that of He II, but there is apparent cospatial 304A and 
I71A emission in the prominence. The explanation for 
this phenomenon is the "prominence-corona transition 
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Fig. 1.— Maps in several SDO/AIA bandpasscs of cavity Dl (a-c) and D2 (f). He II 304A (a). Fe VIII/Fc XXIII 13lA (b), Fc IX/Fc X 
17lA (c), Fe XIV 21lA (d), Fe XII 193A of Dl (e) and D2 (f). Data have radially vignetted to promote structural contrast. 
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region" (PCTR) ( [Orrall fc Schmahl|[l976| . The most 
basic interpretation is tJiat along lines ot sight that in- 
tersect the prominence, there is emission from coro- 
nal, transition region, and condensation plasma. The 
171A emission we see in the prominence is being emitted 
by the transition region between condensed plasma and 
the corona. 

The 193A bandpass shows the cavity very clearly as 
an emission-depletion feature. The coolest mass of the 
prominence can be seen in absorption, which is formed 
by bound-free continuum absorption mostly due to H I. 
For reference, a 193A image of D2 is also presented in 
Figure [if. 

The cavities in Dl and D2 differ primarily in two ways. 
First, D2 is centered around a more angled neutral line. 
It appears to move north to south as the cavity rotates 
onto the solar disk. This has the affect of diffusing the 
contrast between the cavity interior and the streamer. 
Second, the streamer surrounding the cavity in D2 is 
brighter in the coronal bandpasses than the streamer in 
Dl. This can be seen more clearly on the solar disk, 
where the filament channel is significantly more con- 
trasted than in D2. 

2.2. Temporal Statistics on Emission Variability 

Figure [l] served to qualitatively present the morpho- 
logical differences of the prominence cavity in the various 
AIA bandpasses. We will now develop a set of statistics 
to quantitatively describe how the time-dependent cav- 
ity differs from the streamer. The AIA observations can 
be described in the following notation: 

where I is the observed intensity (in data number) in 
bandpass A at time t for CCD pixel ij. We are attempt- 
ing to understand the time dependent plasma processes 
ongoing in the cavity and how those properties differ from 
the streamer, thus we will largely be ignoring the spatial 
coordinates ij. Instead we will subdivide the dataset into 
regions: streamer, cavity, and prominence. The subdi- 
vided regions are shown in Figure [2J Figure [3] illustrates 
how the median intensity (in time at eacn individual 
pixel) in 193A varies between the streamer and the cavity 
as a function of radial height in the plane of the sky. The 
cavity data is presented as red points, while the streamer 
data is presented as black points. The cavity profile is ev- 
erywhere reduced compared to the streamer at the same 
height. This emission depletion varies between 20-50%. 
Assuming an isothermal corona and ignoring projection 
effects, the / n^d£ emission relationship would imply a 
29% densi ty depletion. Th i s is ty pical of cavities as re- 
ported by [Fuller fc Gibson ( 2009 ) . To analyze the time- 
dependent vaHaEiIitynrtEie AIA data, we will use statis- 
tics based on the time series of intensity at every pixel 
within the dataset. 

We have created a statistic, V^^ referred to as variabil- 
ity, designed to quantify the degree to which the intensity 
of a given pixel varies as a function of time. Variability 
is given by 



yij _ ^MAX ^MIN 
^MED 



(1) 



where Z^^/^, ^max' ^^'^ ^med represent the intensity 
value that pixel ij is brighter than through 95%, 5%, 
and 50% of the time series, respectively [I^li^jj will be 
referred to as median intensity). This statistic is akin to 
the variance, however the nature of the intensity distri- 
butions (i.e. asymmetric) does not make variance a par- 
ticularly illustrative statistic. Whereas variance weights 
the entirety of a distribution, our variability statistic 
cares specifically about the "wings" of the intensity his- 
tograms. It measures how the intensity span of the wings 
compares with median intensity. 

Figure [4[ plots variability in the three structural re- 
gions for 17lA, I93A, and 21lA. The upper row of plots 
presents variability versus radial height. Examining the 
streamer distributions for the three bandpasses, we find 
there is a linear relationship between variability and alti- 
tude for all bandpasses. Moreover, the distributions are 
fairly tightly clustered. The cavity and prominence dis- 
tributions are not so simple. The 193A and 211 A data 
both show variabilities which are stronger in the cavity 
than the streamer. The prominence shows higher vari- 
ability than the cavity at low heights. In the 17lA distri- 
butions, there is a striking difference between the cavity 
region and the streamer region. The cavity has a vari- 
ability that is a factor of 4 higher than the variability 
of the streamer. The differences between the variabil- 
ity distributions illustrate that the cavity is much more 
variable than the streamer, but specifically in the cooler 
17lA bandpass. 

The lower row of Figure |4] presents the variability 
data as a function of median intensity. While the up- 
per row illustrated that the cavity is more variable than 
the streamer, it does not describe the significance of that 
variability. To define the significance of the statistic, we 
must define a baseline for what the expected observa- 
tional error should be. With the AIA data, this is chal- 
lenging for two reasons. First, due to the volume of data, 
the AIA datasets are only retrievable as "Level 1" data 
so the engineering level data (flat field and dark current) 
are not available, although these systematics have been 
accounted for at Level to Level 1 conversion. Addi- 
tionally, bandpass imagers cannot photon count because 
there is an indefinite range of photon energies impinging 
on the detector. Due to the non-uniform energy response 
of the detector, this does not allow a singular photon to 
electron to data number gain to be assigned. In the lower 
plots, the expected Poisson noise variability based on the 
preflight detector calibration is plotted as a function of 
intensity as green points. In general, the Poisson vari- 
ability is a close match to the streamer distributions. 
The 211 A data is perhaps the clearest example of this. 
The offset between the streamer distribution and Poisson 
noise can be attributed to a small variability in emission. 

The most important feature to take away from Figurel4] 
is that the variability of the cavity in 17lA is excessively 
strong compared to both the streamer and the baseline 
Poisson noise. There is also excessive variability in the 
cavity compared to the streamer in 193A and 21lA, but 
the I71A data is a factor of 3 more variable. Through 
this statistical analysis, we have identified a unique char- 
acteristic of the EUV emission in the cavity through the 
basis of spatially-independent temporal statistics. Now, 
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Fig. 2. — Time averaged maps in 193A for Dl (left) and D2 (right). The red dashed contour denotes the regions selected as cavity. Blue 
region is prominence. Black region is streamer. These color codes are used in Figures Is] and |4] as well. 
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Fig. 3. — Radial profiles of the median intensity, I*^(p = 0.5), 
for Dl in 193A. As in Figure 2, red points are extracted from the 
cavity region, black points are extracted from the streamer region. 



we will take the next step by honing in on the strongly 
variable regions of 17lA cavity. What is it about the 
intensity distributions in the cavity region that increases 
this statistic? 

Figure [5] shows the histogram for a strongly variable 
cavity point in the 17lA and 193A bandpasses. The vari- 
ability in the 17lA data is over twice that in 193A. In 
terms of the distributions, the most notable feature is the 
strong asymmetry in the 17lA distribution, which has a 
high intensity tail. The 193A histogram, in contrast has 
a symmetric normal distribution. We have chosen statis- 
tics to highlight the nature of these distributions, as the 
intensity histograms plotted in Figurel5]are typical of the 
cavity. The high intensity wing in the 17lA histogram 
make both mean and variance (statistics of specific rele- 
vance in a normal distribution) of little illustrative value. 
Printed along with the variability statistic in Figure [5] is 



Whereas the variability looks at the spread between 
the high-intensity and low-intensity wings of the inten- 
sity histograms, the skew statistic looks at the relative 
intensity-difference between congruent positions in the 
histogram, specifically the brightest 5% of intensities and 
the dimmest 5% of intensities. The skew is a positive 
number, and for the AIA datasets it varies between 0.1 
and 10. A skew of 1 would imply a perfectly symmetric 
distribution. We find the lower cavity is characterized by 
high skew distributions {Bij > 2) in the 171 A bandpass 
but not IQSA bandpass. This signifies that the variabil- 
ity in the 17lA cavity is characterized by short duration, 
strong brightening events. 

Based on the skew and variability statistics, we have 
found that the cavity exhibits very different temporal 
variations in intensity than the streamer. Based on 
these statistics, we now inspect movies of the cavity 
to determine the spatial and temporal coherency of the 
17lA brightening events. 

3. OBSERVATIONS OF PROMINENCE HORNS 

By examining movies of the 17lA bandpass data, the 
nature of these variability features become more appar- 
ent. Figure [6] illustrates a subregion of the cavity going 
through a transient 171 A brightening, and an animated 
version of this figure is presented in the electronic edi- 
tion of this article. The pictured subregion includes the 
prominence and the western half of the cavity. Over the 
course of an hour, we observe the formation of a coUi- 
mated "horn" in the 17lA data which extends nonra- 
dially from the 304A prominence into the lower cavity. 
The apparent width of the feature is 4.2 Mm and the 
apparent length (in the plane of the sky) is 63 Mm. The 
maximum brightening (contrast relative to background) 
is 16%. The spatial geometry of the feature extends from 
the prominence along a concave up loop-segment into the 
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Fig. 4. — Variability statistic plotted against radial height (a-c) and against intensity (d-f) for Dl in 3 different bandpasses: 171A (a and 
d), 193A (b and e), 211A (c and f). The green profiles are the theoretical Poisson-noise variability. Otherwise, the color code is based on 
the regions of Figure 2. 
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Fig. 5. — Intensity histograms for the intensity of a cavity region pixel of Dl in the 171A bandpass (left) and the 193A bandpass (right). 



cavity. 

The correlation between the 17lA feature and the 
structure of the other AIA bandpasses should also be 
considered. In between 2011-02-28 23:24UT and 2011- 
03-01 00:24UT, we find that the prominence emission, 
as seen in the 304A bandpass, has also changed. A co- 
spatial, co-temporal extension of the prominence sits at 
the base of the 17lA horn. The prominence gradually 
extends a maximum of 13 Mm along the 17lA feature 
before gradually contracting. 

This observation demonstrates that there is a coronal 
component associated with prominence dynamics, and 



moreover the dynamic coronal feature projects into the 
cavity interior. These correlated features imply there is 
a magnetic and energetic link between the corona and 
prominence, which leads us to the question: what is the 
physical source for the prominence-corona dynamics? To 
address this question, wc will look towards statistical 
analysis of the many prominence horns to identify phys- 
ical characteristics which will aid us in determining the 
cause and consequence of these features. 



3.1. Method of Analyzing Prominence Horns 
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Fig. 6. — Snapshots at three times in three bandpasses of a prominence horn, zoomed into Dl. Prominence material, He II 304A emission 
(left). Transition region Fe IX-X 171A (middle). Coronal Fe XII 193A (right). The horn is the 17lA brightening which the yellow arrow 
points to. It remains nominally brighter than the background for 1.5 hours. The red contours show the perimeter of strong 304A emission. 
The horn in co-spatial with an extension of the prominence. There is weak, correlated signal in 193A radially above the yellow arrow, but 
overall the structure in 193A is distinct from that of 17lA. A movie of this figure is presented in the electronic version. 
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The first step in this analysis is to build a database of 
events. As we previously established, the most dominant 
forms of emission variability in the cavity are collimated 
regions of strong brightening in 171A emission. We use 
this trait as the defining characteristic of feature identifi- 
cation, and we use visual inspection to identify potential 
features in both Dl and D2. We originally attempted 
to develop an automated method to extracting features 
through a wavelet app roach similar to the OCCULT al- 
gorithm developed by Aschwanden (2010). These tech- 
niques are not well suited to ott-limb corona, where three 
effects hamper them: the strong radial intensity gradi- 
ent, the relatively diffuse structure, and the lower signal- 
to-noise ratio. Instead of investing time into automating 
the initial processing, we have opted to manually identify 
potential features, which are later verified through tests 
described later in this section. We identify a total of 48 
potential features (in 66 total observing hours) between 
our two datasets: 27 features in D2 and 21 features in 
Dl. 

To quantify the changes in intensity as a function of 
time along each feature, we synthesize a curved slit which 
is selected to overlay the potential feature. The data from 
the I93A, I71A, and 304A datasets are extracted in the 
vicinity of the slit and are transformed into a curvilinear 
coordinate system, the axes being distance along the slit, 
y', and distance from slit-center, x' . The transformation 
is done using a triangulation, nearest-neighbor interpo- 
lation scheme. An example of a transformed curvilinear 
dataset is shown in Figure |8] In the temporal domain, we 
extract data from 2 hours prior to the potential feature 
and 3 hours after. 

In order to quantify the time-dependent intensity vari- 
ations in these datasets, we produce light curves as 
a function of position along the slit. The curvilinear 
dataset is divided into background regions (5 pixels< 
\x'\ < 15 pixels) and feature regions (|a;'| < 5 pixels). 
These regions are then coadded at 5 pixel intervals in 
the y'-direction (along the slit). A gaussian temporal 
smoothing is applied to the light curves, with a width 
of 9 minutes. We present light curves which display the 
ratio of feature intensity to the background intensity at 
each segment along the slit as a function of time. Exam- 
ple plots of these light curves are shown in Figure [9| The 
color of each curve denotes in position along the slit. A 
colorbar matching the y' coordinate is plotted in between 
maps in Figure [8] A vertical offset has been applied to 
the light curves shown in Figure [9] so that they do not lie 
on top of each other. 

The final step in the processing of our datasets is find- 
ing and fitting the locations of minima and maxima of 
each individual light curve. This process serves two pur- 
poses: it allows us to automatedly check to make sure 
that the potential features we identified by inspection are 
actually coherent structures, and it allows us to develop 
statistics on the onset, duration, and time-lag between 
different bandpasses of each event. For each light curve, 
we use an algorithm to determine the location of local 
maxima between < i < 40 (where t = 20 is the time of 
maximum feature definition). For a potential feature to 
be labeled as an event, we require that at least 8 sequen- 
tial heights along the slit («17 Mm) have local tempo- 
ral maxima within 3 timesteps of both the neighboring 
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Fig. 7. — A 171A image overlaid with a curved slit. The color of 
the slit represents the curvilinear y-coordinate as shown in Figure 
8. 

heights. This basic test is a simple method which checks 
the coherency of a potential feature. Applying this fil- 
ter to our 48 potential features, we return 45 confirmed 
prominence horns. 

We extend this technique to measure the start and end 
times of each brightening. These time intervals are as- 
sessed by extracting the nearest two local minima on 
either side of the prominence-horn maximum. We apply 
a linear fit to these data (time of event vs. position on 
the slit) where the weight of each minima is inversely 
proportional to its depth relative to the maxima: 

a = {I{t = tuAx) - I{t - tMiN)r^ 

where the linear fit takes the form 

t' — m* y' + tQ. (3) 

The quantity t' is the best-fit time of event (event in 
this context being the initial brightening, peak bright- 
ness, end of subsequent dimming). The derived quantity 
m has units of km~^ s. Its reciprocal is therefore a ve- 
locity which gives us an estimate of how the observed 
disturbance propagates along the slit. The fitting rou- 
tine minimizes the quanity. 



^2^y itobsiy')-t'{y')r 



y' 



(T^ 



(4) 



where tabs is the observed timestep of the event. The 
weighting function values deeper minima more than shal- 
low minima, and strictly-based on number of points, min- 
ima which are well correlated between several heights 
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Fig. 8. — A single horn image, transformed into curvilinear coordinates using a curved slit. In this event, a blob is ejected from the 
prominence to a slit-height of 9 Mm. This ejecta coincides with a coherent 17lA brightening which extends from 5 Mm to 23 Mm. There 
is also a correlated 193A brightening which is not visible in this map because it has not been background-divided. The color table plotted 
horizontally between maps indicate the slit-height of the light curves which are plotted in Figure 9. 



result in lower x^ values. Figure [o] has black lines over- 
plotting the derived linear fits for the feature initiation, 
peak, and dimming. 

3.2. Results 

Using the above described methodology, we will now 
go into characterization of prominence horns. The basic 
configuration for 17lA events is as follows: a rapid, large- 
scale brightening begins with its base near the promi- 
nence. From onset time to peak brightness takes 85 min- 
utes with a standard deviation of 19 minutes. The bright- 
ness enhancement is not constant with distance along the 
feature, and the brightest point along the feature tends to 
occur within the first 15 Mm. The peak brightness con- 
trast generally decreases as a function of distance along 
the horn, such that the horn fades into the background at 
higher heights. We observe that there is a weak correla- 
tion (Pearson coefficient: 0.26) between the peak bright- 
ness of a horn and its length. For the 45 measured horns, 
the median peak brightness contrast is 34%, and 9 fea- 
tures have a peak brightness over 50%. 

Maps showing the positions and strengths of the 
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171A horns are shown in Figure 
trate the proximity of these coronal features to the promi- 
nence. The dense, absorbing regions of the prominence 
appear to be the emanation point. The structure is more 



complex in D2. In the 304A movie, we are able to see 
that there are two unique prominence spines in the D2 
dataset. The smaller northern spine is observed exchang- 
ing plasma with the larger southern spine, which is cen- 
tered under the cavity. These "connection" events have 
similar traits to the more vertical horns in the 17lA light 
curves despite there different geometries. However, they 
do not typically show much structure in 193A. For Dl, 
the overall appearance of the 17lA brightenings are con- 
cave loop segments surrounding a central axis located at 
(+475", +915") at a height of 1.07 R©. 

Of the 45 horns found in the 17lA bandpass, 32 ex- 
hibit correlated brightenings in 193A bandpass data. 
The associated 193A features are plotted in the maps 
of Figure [TT] Of the events which do show correlation 
between 17lA and 193A light curves, there is gener- 
ally not a correlation throughout the entire horn. The 
I93A brightenings tend to occur in the upper sections 
of the I71A horns, the median offset is 6.3 Mm above 
the I71A base. The most striking difference between the 
I71A and I93A features are that the 193A events are 
significantly dimmer. The median peak brightness is 7% 
in I93A compared to the 34% in 17lA. 

From the 32 joint 193A-17lA events, we have 450 co- 
spatial light curves which exhibit correlation. From these 
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Fig. 9. — Light curves of Event 5 (shown in Figure 8) in the 
17lA bandpass (left) and the 193A bandpass (right). The color 
of the light curve denotes its height along the slit (see colorbar, 
Figure 8). The teal curve is extracted from a height of 4 Mm while 
red curve is extracted at a height of 27 Mm. The light curves are 
offset in the vertical direction for clarity. The time of maximum 
brightness contrast of the feature is marked with black crosses, 
while the surrounding minima are marked in black diamonds. The 
linear fit (Equations 3) is overplotted with black lines to show the 
temporal progression of brightening event as a function of slit- 
height. 

points, we measure an average time lag of -0.6 minutes 
±23 minutes between the peak brightness in 17lA and 
the peak brightness in 193A. The same apphes for the 
event onset time. This lack of correlation implies that 
there are both temperature and density perturbations 
causing the brightenings. 

Maps which show the extent of 3G4A features corre- 
lated with 171A horns are shown in Figure 
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„ _,_. We find 
34 of these events. Unlike Figure [lO] and |11[ there is no 
color range applied to the 3G4A features. We did not 
extend the same intensity analysis to the 304A dataset 
because of the large intensity gradients on small-spatial 
scales within the prominence. The basis of our compar- 
ison with 304A data is changes in the perimeter of the 
prominence, as was seen in Figure |6] There are generally 
two type of events we observe in the 304A data. In 20 
of the 304A events, we observe simple extensions of the 
prominence. In the other 14 events, we actually observe 
ejection of prominence material such that no 304A emis- 
sion is detected between the 304A ejecta and the promi- 
nence itself. In both cases, the median rise time and fall 
time is 24 minutes. The median maximum length for 
ejecta is 10.5 Mm, while the median maximum length 
for extensions is 7.6 Mm. In general the rise of 304A ma- 
terial begins after the initiation of the 17lA brightening 
by 30 minutes. We find a Pearson coefficient of 0.44 for 
the correlation of the peak 304A extension time and the 
time of peak 17lA brightness. 



Velocities are a fundamental quantity associated with 
dynamics, however this is a difficult measurement to 
make in prominence horns. We do not have spectral 
Doppler observations for these cavities, and the applica- 
tion of feature tracking techniques is not obvious; our fea- 
tures appear to maintain a constant spatial orientation 
and position over their 3 hour duration. As described in 
Section 3, we have attempted to assess a velocity based 
on the time of initial brightening and peak brightening. 
While this quantity has physical units, we must be clear 
on what is being measured in terms of basic physics: 
there is a perturbation happening along a coronal loop 
which is changing the amount of photons being emitted 
along that loop as a function of tim,e. We are measuring 
the time that a discernible change in photoemission has 
occurred at a particular position along the loop relative 
to other positions on the loop. 

The derived velocities based on the initial brightening 
and the peak brightening are shown in Figure [13] as a 
function of feature length and reduced x^ of the model 
fit (using Equations 3 and 4). Feature length is impor- 
tant in this diagnostic because we have more points to 
track the apparent propagation of the feature. The pe- 
culiar X-axis scale is due to the nature of the fitting. The 
fit-parameter, m is converted into a velocity by 

Ay' _i 

V = — ; — m 
At 

where Ay' represents the plane of sky distance of 5 slit- 
pixels (2100 km) and At is the cadence of the data (360 
s). A negative value of m would imply a propagation 
from higher radial heights downward, while a positive m 
implies a propagation from the prominence upward. The 
median derived m for the 17lA data for the brightness 
initiation is 0.056 which corresponds to a velocity of 104 
km s~^. Given the errors in measurement, we are not 
able to discriminate between high speed (|w| > 100 km 
s~^) features moving upward and downward. 

The derived velocities in Figure [13] are sparsely dis- 
tributed, but the highest quality fits tend to cluster be- 
tween \v\ > 20 km s""'^. These velocities are intriguing 
considering the event durations. A transient disturbance 
traveling at 20 km s~^ would advance 60 Mm in 5 min- 
utes. Our observed features last longer than 1 hour. This 
implies that the propagation of the perturbation and the 
subsequent equilibration or damping of the perturbation 
are based on different energetic processes and timescales. 

3.3. Estimate of Horn Density 

We wish to estimate the density of the horns, and this 
re quires us to account for p roje ction effects. Past e fforts 



by [Schmit fc Gibson' (2011) and ]Kuceraet al.] (l2012') have 
gone through modeling ettbrts to quantify the density and 
temperature effects which are manifest in optically-thin 
EUV emission. Extending this work to the Dl cavity, we 
can make a back-of-the-envelope calculation for promi- 
nence horns. Let us consider a line of sight which ob- 
serves a 53° long cavity embedded with in a 70° long 
streamer at a latitude of 60°. Using the forward mod- 



eling technique of Gibson et al. (2010), we predict there 
should be approximately 9 x 10^*^ cm~^ of emission mea- 
sure along the cavity-center line of sight at a height of 
1.07Rs. Assuming an isothermal cavity and streamer of 
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Fig. 10. — 193a map of Dl (left) and D2 (right) in grayscale. Location of horns in the 171A bandpass. Color table advances black to 
orange to white with increased maximum (in time) contrast to background. Range: 4% to 50% 
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Fig. 11. — 193A map of Dl (left) and D2 (right) in grayscale. Location of brightenings in the 193A bandpass associated with horns. 
Color table advances black to blue to white with increased maximum (in time) contrast to background. Range: 1% to 20%. 



1.7 MK, we can write the intensity ratio of the promi- 
nence horn to the static cavity in the form, 

I'm ^ MiijT = 1.7MK) ^ EM, +^in(r = Tpeafc) * EM^ 
/i7i Am{T =l.7MK)*¥M 

(5) 
where A represents the AIA temperature response func- 
tion per unit emission measure. The forward model pre- 
dicts that most of the ambient 17lA emission comes from 



low-density, weakly-emitting plasma, but the plasma is 
quite extended (the full-width half-max for emissivity 
distribution is ±70 Mm out of the plane of the sky). 
The ratio of Ani{T = 1.7 MK)/Ani{T = Tp^ak) is ap- 
proximately 4.5%. We have divided emission measure 
distribution into three components 



EM = / n^ dl, 
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Fig. 12. — 193A map of Dl (left) and D2 (right) in grayscale. Location of prominence extensions in 304A emission associated with the 
horns. These tracks represent the range of prominence extension or ejection. No color range is applied. 



Event Start 



100 - 



80 - 



E 

^ 60 

c 

I 40 

D 
<U 

L- 

20 



+ + : ++ 
+ 












Event Pea 


k 










100 


- 














- 


80 


- 




x + 










- 


60 


~ 














~ 


40 


_ 




+ + Hi * 










_ 


20 


- 


+ 
+ 


+ : + + + 
+ : + 
++ + + 

+ : 


t 


+ 
+ 


+ 


+ 
+ 





-10 



-25 -100 100 25 10 

Projected Velocity (km s'') 



05 



-10 



-25 -100 100 25 10 

Projected Velocity (km s"') 



05 



Fig. 13. — Scatter plots of the derived horn velocity along the slit versus horn length. Color table advances black to orange to white 
with increasing x^ value. Derived velocity for the horn initiation (left) and horn maximum (right). Black crosses represent the most 
reliable velocities, x^ range for initiation: 0.1 to 1.4; x^ range for maximum: 0.8 to 6.5. The black dotted line represents an infinite and 
directionless speed (m = 0). 



to produce a 50% increase in 17lA intensity 



EM., = ; n- 

'3 Mm 



n{i = 0) 



1.2 



EM' 



3 Mr 



where n' is the horn density and n is radial density pro- 
files from Schmit & Gibson (2011). These emission mea- 
sures are the full line of siglit, the horn region in the am- 
bient cavity, and the horn region during the 171A bright- 
ening respectively. We estimate the line-of-sight width of 
the horn at 6 Mm. Substituting in these equations, we 
can solve the minimum density perturbation necessary 



This increase in density assumes that the plasma has also 
cooled to 1 MK, where it emits 18 times more strongly 
than at 1.7 MK, the original imposed temperature. This 
assumption makes it such that the predicted density 
change is a minimum value, whereas a weaker tempera- 
ture perturbation would result in a weaker enhancement 
in emission. Overall, we predict that a 6 Mm wide swath 
of plasma would have to increase its emission by a factor 
of 22 to produce an overall brightening of 50% when ac- 
counting for background emission. 

We can apply this same technique for 193A emission 
where we use the same temperatures and the derived n' 
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and solve for 1(93 //iga. We now have the ratio in the 
temperature response function of 



AMT = 1.7 MK) 
AMT = 1 MK) 



= 2.4 



The reduction in temperature which raised the emissivity 
for the I71A bandpass by a factor of 20, acts to reduce 
the I93A emissivity. The predicted ratio I'lg^/Iigs has a 
value 0.98 compared to 1.5 for I[j^/Ini. In other words, 
the intensity produced by a corona dominated by > 1.5 
MK structure will not vary drastically for localized (short 
line of sight contributions) cool-coronal events. 

4. CONCLUSIONS 

We have demonstrated that the cavity exhibits fun- 
damentally stronger dynamics than the surrounding 
streamer. This was illustrated through a statistical anal- 
ysis of EUV observations from the SDO/AIA instru- 
ment. In particular, the statistics indicate that the 
lower-cavity region surrounding the prominence under- 
goes strong brightening events in the 17lA bandpass, 
which is typically dominated by 1 MK plasma. 

Analyzing movies in the 17lA bandpass, the obser- 
vations indicate that these events are coronal loop seg- 
ments, called prominence horns, which extend from the 
prominence into the corona. Horns project into the 
interior of the cavity. The 17lA emission of horns 
is often correlated with weakly-enhanced emission in 
I93A bandpass, which contains hotter coronal emission. 
The I71A emission is often co-spatial with extensions of 
the prominence material as seen in the 304 A bandpass. 

We have gone through many details of these emission 
signatures, but we have not yet addressed the physical 
nature of these events. There are several obvious puz- 
zles present in the data. First, how do we explain the 
apparent co-spatial 304A and 171 A emission? There is 
a large temperature discrepancy between the theoreti- 
cal peak emissivity of He II 304. 3A and Fe IX 171. 4A. 
There are two possible explanations: cool blends in the 
I71A bandpass or projection effects. 

There are cool blends in the red periphery of the 
I71A bandpass, in parti cular resonance li n es of O V 
I72.2A and O IVI73.OA. 
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(2011) con 



eluded thorough analysis ot Hinode/t)lS data that the 
contributions of these lines is limited to a few percent on 
disk. These lines will not emit strongly in equilibrium 
because of the large difference between the excitation 
energy (70 eV) and the temperature associated with the 
peak ion population (O VI population peaks at 4x10^ 
K « 35 eV). Radiative recombination could play a role 
in populating these excited levels in a non-equilibrium 
processes. However, the magnitude and duration of this 
recombination emission is unlikely to explain prominence 
horns. 

Another possible explanation for the correlated 
3O4A and I71A emission is projection effects. Horns may 
have a significant component of their geometry aligned 
along the line of sight. Emission which appears cospa- 
tial would occur along separate segments of an individual 
magnetic field line. The prominence-corona transition re- 
gion, which would separate these segments, is predicted 
to be less than 5 Mm thick. 



Another peculiar aspect of prominence horns is the 
weak correlation between 17lA emission and 193A emis- 
sion. There are relative weak separations between the 
ionization energies for Fe IX and Fe XII. This has 
the effect that the predicted AIA temperature response 
function has significant overlap between the 17lA and 
193 A bandpasses. Moreover, there is precedence that Fe 
IX 171 A and Fe XII 195A emission are often highly corre- 



lated in dynamics heating events (jViall fc Klimchuk 2011 



Aschwanden et al. 2000). The most likely explanation 



tor the weak correlation in horns is that the dominant 
component emission in the 17lA is formed at tempera- 
tures less than 0.8 MK. This cool plasma would not have 
strong emission in the 193A bandpass. 

In Section 2, we discussed that there was more variabil- 
ity in the 17lA bandpass than the 193A bandpass. While 
horns explain our statistical analysis of the 171 A cavity, 
the associated 193A component of the horns does not ex- 
plain the variability in that bandpass. In examining the 
light curves, we find that there are intensity variations 
a factor of 3 larger in the 193A cavity than the inten- 
sity variations related to the horns. The power seems to 
be strongest for variations with a 4-8 hour period and is 
more evenly distributed between dimming and brighten- 
ing event s. We believe that it is t his variability compo- 
nent that Wang & Stenborg (2010) have studied through 
their wavelet analysis. 

We have created a framework with which we could 
measure the characteristic velocity of the prominence 
horns, but found that the spread in the derived values 
was large relative to the values themselves. We find 
several examples from our 45 events where an entire 50 
Mm long structure undergoes a uniform and simultane- 
ous change in intensity. In terms of basic physics, these 
large-scale isotropic changes are most readily explained 
through thermal conduction. 

Conduction is a diffusive operator. Given a thermal 
perturbation on a coronal loop, anisotropic thermal con- 
duction rapidly redistributes the energy associated with 
the perturbation. These characteristics match the scale 
(> 50 Mm) and velocity (> 50 km s~^) of horns. In the 
corona however, thermal changes and density changes 
are inseparable. Due to the different velocity scales as- 
sociated with thermal conduction and siphon flows, a 
change in temperature will always result in a change in 
density. A change in density, given mass conservation, 
requires velocities. The hydrodynamic equations must 
be self-consistently solved to truly understand the impli- 
cations for any thermodynamic perturbation. Of special 
interest to these observations is a detailed analysis of a 
catastrophic cooling model for prominence formation. A 
followup paper is underway comparing the light curves 
of promi nence horns with th e thermal non-equilibrium 
model of Karpen et al. ( 2006 ) . 

We have not yet discussed the observed geometry 
of the prominence horns. Horns are tightly-collimated 
and exhibit little structural variation over their du- 
ration. These characteristics indicate that these are 
magnetically-aligned structures, or more specifically seg- 
ments of coronal loops. We have speculated they may be 
explained as thermal structures, and thermal conduction 
is dominantly field aligned. Thus, we posit the extracted 
prominence horns are indicative of the magnetic field ge- 
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onietry of the lower cavity. 

In both Dl and D2, a discernible axis can be found 
which the horns circumscribe. Horns imply that there 
is a component of the prominence magnetic field which 
is non-axial. This is an indication that twist is present 
in the quiescent prominence-cavity system. However, we 
cannot be completely certain that these magnetic struc- 
tures are within the cavity itself. It remains possible that 
horns and the cavity inhabit separate field lines which 
overlay each other in optically thin projection. It is also 
interesting that despite the relative symmetry of horns 
about the prominence spine, we do not find evidence that 
horns on one side are correlated with horns on the other 
(we do not observe dynamics which extend through the 
prominence spine). 



These observations are the first attempt to statisti- 
cally quantify the connection between cavity dynamics 
and prominence dynamics. We have identified a strong 
emission variability connection between the prominence 
and the cavity, but we have not gone so far as to deter- 
mine the physical source of these dynamics. There is a 
compelling paradigm for the prominence-cavity connec- 
tion: a mass exchange feeds condensing plasma out of the 
low-density cavity into the cool, overdense prominence. 
The models for catastrophic cooling make predictions on 
what coronal perturbations are necessary to drive this 
process. In the forthcoming article, we will test these 
predictions against the emission variability we have doc- 
umented in the prominence environment. 
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